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The  optical  properties  of  particles  deposited 
on  a  surface. 

Final  Teclinical  Report  on  Contract  DAJA45-93-C-0043 

1.  Outline  of  the  research. 

In  recent  years  the  optical  properties  of  particles  near  or  on  a  plane  surface 
have  been  the  subject  of  several  papers  that  deal  with  both  the  theoretical  and 
the  experimental  aspects  of  the  problem.  This  interest  from  the  scientific  com¬ 
munity  can  be  easily  understood  because  revealing  through  optical  measure¬ 
ments  the  presence  of  particles  on  a  surface  that  is  expected  to  be  clean  is  of 
great  importance  both  for  industry  and  fc  nedicine. 

Among  the  theoretical  techniques  thr  “  been  proposed  to  deal  with  the 
problem  the  method  of  the  images  prov  ^  the  most  fruitful  when  the 

plane  surface  is  a  perfectly  reflecting  one,  a  .i^’■.)n  that  is  met  to  a  high  de¬ 
gree  of  accuracy  by  many  polished  metallic  sur;  .s.  According  to  the  theory  of 
the  images  the  problem  of  scattering  from  a  particle  in  the  vicinity  of  a  perfectly 
reflecting  plane  surface  is  quite  equivalent  to  the  problem  of  the  dependent 
scattering  from  the  compound  object  that  is  composed  by  the  actual  particle 
and  by  its  mirror  image  when  excited  by  the  superposition  of  the  fields  that 
come  from  the  actual  source  and  from  its  image.  Of  course,  the  method  of  the 
images  is  advantageous  provided  one  is  able  to  solve,  without  undue  computa¬ 
tional  effort,  the  problem  of  scattering  from  the  compound  scatterer.  In  fact,  the 
lack  of  an  efficient  technique  to  deal  with  this  problem  forced  several  research¬ 
ers  to  solve  the  problem  for  limiting  cases  only,  e.  g.  normal  incidence,  or 
through  approximations,  such  as  the  assumption  of  smallness  of  the  scattering 
particle. 

2.  Method  of  attack. 

The  method  that  we  used  for  our  research  combines  the  method  of  the  imag¬ 
es  with  the  technique  that  we  developed  several  years  ago  to  calculate  the  scat¬ 
tering  properties  of  aggregated  spheres.  When  one  assumes  that  the  actual  par¬ 
ticle  is  or  can  be  modelled  as  a  cluster  of  spherical  scatterers,  the  method  of  the 
images  requires  to  deal  with  the  properties  of  a  cluster  that  contains  twice  as 
many  spheres  as  the  actual  object.  The  fact  that  the  exciting  field  is  the  superpo¬ 
sition  of  two  plane  waves  does  not  introduce  any  additional  difficulty.  Further¬ 
more,  the  presence  of  the  reflecting  surface  does  not  prevent  us  from  perform¬ 
ing  analytical  averages  over  the  orientations  of  the  scattering  particles. 
Therefore  our  approach  is  suitable  to  calculate  the  scattering  pattern  not  only 
from  single  spheres  and  from  single  aggregates  of  spheres  but  also  from  a  dis¬ 
persion  of  identical  aggregates  whose  orientations  differ  only  for  a  rotation 
around  an  axis  orthogonal  to  the  reflecting  surface. 

As  our  approach  does  not  imply  any  approximation  nor  limitation  both  on 
the  direction  of  incidence  and  of  observation  and  on  the  polarization  both  of  the 
incident  and  of  the  scattered  field,  we  were  able  to  calculate  the  full  scattering 
pattern  for  a  few  representative  examples  of  aggregates  that  are  meant  to  simu- 


late  the  presence  of  anisotropic  particles  on  the  reflecting  surface.  Moreover, 
our  approach  applies  also  to  hemispheres  and  to  aggregated  hemispheres  with 
their  flat  face  lying  on  the  surface;  this  suggests  the  possibility  of  simulating  the 
properties  of  an  arbitrarily  rough  metallic  surface. 

3.  Sketch  of  the  results. 

Although  both  the  theory  and  the  main  results  of  our  research  are  fully  ex¬ 
plained  in  the  enclosed  paper,  that  has  been  accepted  for  publication  in  JOSA 
A,  it  may  be  useful  to  summarize  here  the  most  relevant  findings. 

For  a  dispersion  of  model  anisotropic  particles  on  the  surface  our  results 
show  that  the  anisotropy  is  not  masked  by  the  random  orientation  of  the  parti¬ 
cles;  the  anisotropy  can  be  revealed  by  observing  the  polarized  light  that  propa¬ 
gates  along  the  reflecting  surface  at  right  angles  to  the  plane  of  incidence. 

The  need  to  deal  with  the  compound  scatterer  (actual  particle  and  its  image) 
persists  until  the  distance  between  the  actual  particle  and  the  surface  does  not 
exceed  5d,  where  d  is  the  linear  size  of  the  actual  particle  (for  spheres,  d  coin¬ 
cides  with  the  diameter). 

The  approximations  that  are  based  on  the  assumption  of  smallness  of  the 
scattering  particles  may  present  severe  drawbacks.  In  fact,  although  a  small 
particle  in  the  vicinity  of  the  surface  is  equivalent  to  the  cluster  of  two  identical 
small  particles,  the  cluster  itself  may  be  not  small  enough  to  grant  the  applica¬ 
bility  of  the  Rayleigh  approximation. 
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In  the  framework  of  the  image  theory  (he  full  scaltcriiig  pattern  from  model  anisotropic  particles 
on  a  perfectly  reflecting  surface  is  calculated  for  an  arbitrary  direclion  of  the  incident  field.  The 
particles  on  the  surface  are  modelled  either  as  clusters  of  spherical  scatterers  or  as  clusters  of  hemi¬ 
spheres  who,se  flat  face  lies  on  the  reflecling  surface.  Oiir  approach  is  ba-sed  on  the  e.vpansion  of  all 
the  fields  in  terms  of  spherical  multipolcs  whose  transformation  properties  are  used  to  get  a  compact 
expression  for  the  scattered  intensity  both  from  a  single  particle  as  a  function  of  its  orientat'on  and 
from  a  dispersion  of  randomly  oriented  particles.  The  patterns  that  were  calculated  for  several  model 
scatterers  show  some  features  that  may  give  useful  suggestions  on  the  possible  anisotropy  of  actual 
particles  on  a  reflecting  surface. 


1.  Introduction 

The  exact  solution  to  the  problem  of  electromagnetic  scattering  from  particles  near  or  on  a  substrate  of  general 
dielectric  properties  bounded  by  a  plane  surface  has  been  worked  out  for  spheres  only.'"'*  The  literature  also  reports 
approximate  solutions  both  for  spheres  and  for  particles  of  general  shape  that  apply  e.  g.  to  objects  far  from  the 
surface  or  to  small  objects  on  the  surface  itself.®''^  The  need  to  resort  to  aj>proximations  is  not  surprising  because  the 
simultaneous  presence  of  the  closed  surface  of  the  particle  and  of  the  open  plane  surface  of  the  substrate  complicates 
the  boundary  condition  problem  that  one  has  to  solve  to  get  the  pattern  of  the  scattered  radiation.  When  the  plane 
surface  is  perfectly  reflecting  the  image  theory  allows  one  to  substitute  the  problem  at  hand  with  the  equivalent 
problem  of  scattering  from  the  object  that  is  composed  by  the  actual  .scatterer  and  by  its  image/  but  this  substi¬ 
tution  amounts  to  an  effective  simplification  only  when  the  dependent  scattering  from  the  compound  object  can  be 
calculated.  Actually,  in  the  case  of  a  perfectly  reflecting  surface  an  exact  solution  has  been  worked  out  for  both 
spheres®’*  and  cylinders.® 

In  this  paper  we  will  resort  just  to  the  image  theory  to  calculate  the  pattern  of  the  scattered  radiation  from  single 
model  anisotropic  particles  on  or  near  a  perfectly  reflecting  plane  as  well  as  the  pattern  from  a  dispersion  of  identical 
model  anisotropic  particles,  with  random  orientat  ion,  that  form  a  low-density  monolayer  on  or  near  the  surface.  The 
anisotropic  particles  will  be  modelled  as  clusters  of  spherical  scatterers  on  account  that,  in  the  framework  of  the 
image  theory,  the  scattering  from  a  single  sphere  near  a  perfectly  reflecting  surface  is  equivalent  to  the  dependent 
scattering  from  the  binary  cluster  that  is  composed  of  the  actual  sphere  and  of  its  image,  provided  the  exciting  field 
is  the  superposition  of  the  fields  that  come  from  the  actual  source  and  from  its  image.*  Thus,  if  the  object  in  front 
of  the  surface  is  a  cluster  of  spheres,  one  has  to  deal  with  the  .scattering  pattern  from  a  cluster  that  has  twice  the 
spheres  than  the  actual  cluster. 

The  scattering  from  a  cluster  of  .s|)heres  ran  he  calculated  with  a  high  degree  of  precision  and  without  undue 
computational  effort  through  the  tcchnifiue  that  wc  devised  a  few  years  ago.'®  We  based  our  approach  on  the 
expansion  of  the  incident,  of  the  scatlered  and  of  I  he  internal  field  in  terms  of  spherical  multipoles"”''^  and  on  the 
imposition  of  the  appropriate  boundary  coii<lit ions  at  the  surface  of  each  cnie  of  the  spheres  in  the  cluster.  The 


final  equations  are  afiected  by  no  approximation  and,  in  fact,  prove<l  to  yield  results  that  are  in  excellent  agreeirient 
with  the  available  experimental  data  for  single  clusters.'®  '®  A  noticeable  feature  of  our  formalism  is  that,  through 
the  use  of  the  transformation  properties  of  the  spherical  multipoh's  under  rotation,  the  scattering  amplitude  of  any 
cluster  turns  out  to  be  factori/ed  into  a  part  that  depends  only  on  the  structure  and  on  the  orientation  of  the  cluster 
and  into  another  part  that  depends  only  on  the  polarisation  and  on  the  direction  both  of  the  incident  and  of  the 
scattered  field. On  account  of  this  separation  the  pattern  of  the  scattered  radiation  from  a  cluster  of  arbitrary 
orientation  can  be  calculated,  without  undue  computational  clfort,  for  whatever  polarization  and  direction  both  of 
incidence  and  of  observation.  Our  formalism  is  also  suitable  to  calculate  the  scattering  pattern  from  an  assembly 
of  identical  clusters  with  random  orientation,  provided  their  number  density  is  so  low  that  the  multiple  .scattering 
processes  among  different  clusters  can  be  neglected.  The  scattering  amplitude  depends,  indeed,  on  the  orientation 
of  the  cluster  only  through  the  D-matrices,  i.  e.  through  the  irreducible  representations  of  the  three-dimensional 
rotation  group;'®  therefore,  once  the  scattering  amplitude  is  known  for  a  single  orientation,  it  is  al.so  immediately 
known  for  any  other  orientation.  Using  this  dependence  of  the  scattering  amplitude  on  the  orientation,  the  average 
that  must  be  calculated  to  get  the  scattering  pattern  from  the  whole  as.sembly  can  be  performed  analytically  in  many 
cases  of  interest;  even  when  the  average  must  be  performerl  numerically  the  computation  requires,  however,  little 
effort. 

In  this  paper  we  will  use  ju.si  the  formali.sm  that  we  outlined  above"’-' ‘  to  calculate  the  .scattered  intensity  from  a 
dispersion  of  identical  clusters  on  the  reflecting  surface;  our  only  assumption  is  that  the  actual  objects  together  with 
their  images  are  still  identical  to  each  other  and  that  the  orientations  of  any  two  clusters  differ  for  a  rotation  around 
an  axis  orthogonal  to  the  refleci  iug  surface.  We  stri*ss  that  our  procedure  does  not  imply  any  limitation  either  on  the 
polarization  or  on  the  direction  both  of  the  incident  and  of  the  observed  wave.  The  same  procedure  is  also  applicable 
to  one  or  naore  hemispheres  whose  flat  face  lies  on  the  surface  as  the.se  objects,  together  with  their  images,  form 
either  a  single  sphere  or  a  cluster  of  spheres.  In  this  paper  wc  will  also  deal  with  the  scattering  pattern  from  these 
objects  in  order  to  gain  some  information  on  their  properties  because,  in  our  opinion,  a  dispersion  of  hemispheres  on 
a  plane  surface  could  be  an  useful  motlel  to  simulate  the  properties  of  e.  g.  a  rough  metallic  surface. 


2.  Theory 

A.  The  incident  ficbl  and  the  i-eflocted  flehl 

Let  us  assume  that  the  halfspace  ;  <  0  is  filled  with  a  perfectly  conducting  material,  such  a.s  a  metal:  the  bounding 
surface  z  =  0  is  thus  perfectly  reflecting.  We  a.ssume  that  the  incident  field  is  the  plane  wave 

E’"'  =  Eoc,  exp(?k,  r) 

that  propagates  in  the  halfspace  z  >  0  of  (real)  refractive  index  n  with  (unit)  polarization  vector  e/  and  wavevector 
k/  =  nJbic/;  as  usual,  k  =  w/c  and  the  time  dependence  exp(— «*;<)  is  omitted  throughout.  The  multipole  expansion 
of  the  incident  field  is 

plm 

where  we  define  the  spherical  multipoles 

=  i/(itr)X,„,(f),  X  iyj(k,v) 

whose  parity  is  given  by  the  superscript  p  =  1,2  that  distinguishes  the  magnetic  multipoles  (p  =  1)  from  the  electric 
ones  (p  =  2);  the  X;,,,  s  are  vector  spherical  harmonics.'''  In  turn  the  mnltipole  amplitudes  Wj  are  defined  as 


where 


k)  =  Ijr.'c  x;„.(k),  \vll\e,  k)  =  (k  x  c)  X;,„(k). 

The  reflected  field  coincides  with  the  field  that  conies  from  the  mirror  image  of  the  source  and  has  the  expansion 

plm 

Of  course,  the  condition  of  reflection  implies  that  the  polarizations  of  E'"'  and  E"-^  cannot  be  mutually  independent; 
therefore  the  relation  between  them  will  be  presently  established.  Since  the  plane  z  =  0  is  perfectly  reflecting,  there 
is  no  transmitted  field:  thus  E"*'  and  E"-^  must  have  the  same  amplitude  as  well  as  the  same  phase  on  that  plane. 
Therefore,  the  boundary  conditions  at  the  reflecting  surface  read 

(e/  +  e«)  X  1*1  =  0,  (2) 

where  n  is  the  unit  normal  to  the  surface.  In  order  to  find  the  explicit  relation  between  the  polarization  of  the 
incident  and  of  the  reflected  wave  it  is  convenient  to  introduce  a  basis  for  the  projection  of  the  states  of  polarization. 
For  the  incident  wave  we  introduce  the  pair  of  unit  vectors  ii/||,  that  lies  in  the  plane  of  k/  and  of  fi  and  is  orthogonal 
to  k/,  and  u/j.  =  k/  x  u/||;  for  the  reflected  wave  we  define  the  pair  u/}||  and  u/jx  that  are  orthogonal  to  k/t  and  to 
each  other  with  uh||  in  the  plane  of  k/j  and  u.  By  projecting  6/  and  cn  on  the  respective  basis  and  imposing  the 
boundary  conditions,  e<i,  (2),  we  get 

eftll  =  e/||.  e/ti  = 

when  the  polarization  is  linear.  For  the  ca.se  of  circular  polarization  we  define  the  ba.scs 
T  'ti/x).  »/?T  = 

and  with  the  same  procedure  we  get 

=  «/+.  =  e/.. 

Although  the  field  that  actually  excites  the  particles  is  E'”  =  E’"'  +  E’*-^,  the  relations  above  allow  us  to  refer  to 
the  direction  of  incidence  and  to  the  polarization  of  the  incident  field  only. 


(uflil  qp 


D.  The  scattered  field 

The  field  that  is  scattered  by  the  compound  object  that  includes  the  actual  particle  and  its  mirror  image  can  be 
expanded  as 

E*«  =  (3) 

plm 

where  the  multipole  fields  HJ,^,*’s  are  identical  to  the  jj,'^’’s  except  for  the  substitution  of  the  spherical  Hankel  functions 
h\^\kr)  to  the  spherical  Be-s-scl  functions  ji(lcT).  The  amplitudes  are  related  to  the  amplitudes  of  the  exciting 
field  E“'  through  the  equation 

a(p)  _  _  V'  (4) 

•^/m  ~  El'm” 

p'/'m' 


ixAp')  _  ii/lf''  j. 

'''  EVm'  —  H'm'  ^  ’^ni'm" 
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The  matrix  S,  of  eloim  iils  is  the  so  calk'd  transition  matrix'^  and  contains  all  the  stnictnral  information 

on  the  aggregate  composed  of  the  particle  and  of  its  image.  The  form  of  ecp  (i)  is  essential  to  perform  the  angular 
analysis  of  the  scattered  radiation,  as  in  this  etpialion  the  slriictiiral  information  is  factorized  with  respect  to  the 
angular  information.  In  other  words,  S  depends  on  the  geometry,  on  the  orientation  and  on  the  scattering  power  of 
the  particle  whereas  all  I  he  information  on  the  direction  ainl  on  the  polarization  of  the  incident  wave  is  contained  in 
the  We’s. 

Let  us  now  recall  thai  the  particles  that  we  deal  with  in  this  paper  together  with  their  image  always  form  a  cluster 
of  spheres;  therefore  their  S-matrix  can  he  evaluated  through  the  technique  that  we  devised  several  years  ago  to 
describe  the  scattering  properties  of  aggregated  spheres.  The  complete  theory  has  been  |)ublished  elsewhere*®  so 
that,  in  Appendix  A,  we  report  only  the  guidelines  of  our  procedure  for  the  calculation  of  S.  This  matrix  is  related 
to  the  matrix  .M  *  that  is  the  inverse  of  the  matrix  .'td  that  e.xpresst's  the  boundary  conditions  at  the  surface  of 
each  one  of  the  spheres  in  terms  of  multipole  fields.*'*  The  matrix  M  accounts  for  the  fact  that  the  incident  field 
on  each  one  of  the  spheres  is  the  superposition  of  and  of  the  fields  that  have  already  been  scattered  by  all 
the  other  spheres  in  the  aggregate.  In  other  words  the  field  scattered  by  the  whole  aggregate  is  calculated  through 
a  dependent  scattering  approach.  For  an  aggregate  of  N  spheres  with  no  particular  symmetry  the  order  of  M  is 
Ja/  =  2N Lm(Lsi  +2),  where  La/  is  the  maximum  value  of  /  that  must  be  retained  into  the  expansions  eq.  (1)  and  (3) 
to  get  fairly  convergent  results  for  the  scattered  field.  Therefore,  .M  can  become  rather  big  thus  requiring  a  long 
computer  time  for  its  inversion;  actually,  the  calcidation  of  A-f'*  is  responsible  for  the  most  part  of  the  calculations 
that  are  needed  to  get  the  sc, altered  field.  The  us<'  of  group  theory,  however,  may  help  to  reduce  the  computational 
effort  through  the  exploitat  ion  of  tin?  .symmetry  |)roperties  of  the  scatterer."®'"*  In  this  respect,  it  is  worth  recalling 
that  our  compound  .sratterers  (actual  object  and  its  image)  do  have  at  least  the  reflection  symmetry  with  respect  to 
the  plane  of  the  substrate.  Furth('rmore,  if  the  compound  .scatterers  have  also  the  cylindrical  .symmetry  around  an 
axis  orthogonal  to  the  snrfar(‘.  the  dio'ct  use  of  the  machinery  of  group  theory  is  not  necessary,  because  the  choice 
of  the  cylindrical  axis  as  the  z  axis  produces  an  automatic  factorization  of  the  .Ad  matrix  with  respect  to  the  index 
m.  In  this  ca.se  the  highest  oiah'r  of  the  tnatrices  to  be  inverted  is  dst  —  2N  L,\i  only. 


C.  The  scattered  iiiteii.sity 


The  relevant  quantity  is  the  matrix  /,,,,<  =  ;  the  indices  and  /?'  refer  to  the  polarization  of  the  scattered 

and  of  the  incident  field.  r('specl,ively.  'I'he  scattered  field  from  a  particle  at  the  origin  of  the  frame  of  reference  can 
be  written  as 


r5ca  _  r'  exp(l«tf)  , 


provided  r  is  large,  so  that 

l2 


fl)!)'  — 


/t')' 


where  the  /„»  are  the  elements  of  the  normalized  scattering  amplitude*"*  that  are  known  to  depend  both  on  the 
incident  and  on  the  scattered  wavevector  as  well  as  on  the  orientation  of  the  particle  with  respect  to  a  frame  of 
reference  fixed  in  the  laboratory.  Of  course,  the  scattering  amplitude  that  we  must  consider  here  is  that  of  the 
compound  scatterer  that  is  composed  by  the  object  in  front  of  the  surface  and  by  its  mirror  image.  The  /,,,<  are 
related  to  the  elements  of  the  transition  matrix  through  the  equation*" 


Inn'  — 


Aitnk 


plm 


’  *  ,s' 


•^hn,rw'  ICri't'w' 


(5) 


with 


where  es  and  ks  are  the  polarization  vector  and  the  direction  of  observation  of  the  .scattered  field,  respectively.  If 
one  has  to  calculate  the  scattering  pattern  from  a  single  cln.ster  in  the  vicinity  of  the  reflecting  surface,  eq.  (.5)  is 


quite  efficient  as  i>s  form  shows  that  ,  onrc  i  lic  S-iiia(rix  is  known,  the  scattered  intensity  can  be  calculated  without 
undue  computational  effort  for  any  k/  and  k^. 

When  one  is  interested  in  the  pattern  of  t  he  scattered  radiation  from  an  assembly  of  identical  scatterers  in  front 
of  the  surface,  the  distribution  of  their  orientations  must  be  taken  into  account.  This  can  be  done  by  associating 
to  each  cluster  a  local  system  of  axes  chosen  so  that  when  two  clusters  are  stiperposetl  their  respective  local  axes 
coincide.  As  explained  el.s<.'where  in  full  detail, the  scattering  amplitude  of  a  cluster  of  known  orientation  is  given 
by 
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where  ©  =  (0',/j, 7)  is  a  shorthand  for  the  three  Euler  angles'"  that  characterize  the  orientation  of  the  local  frame 
with  respect  to  the  laboratory  frame.  In  etp  (ti)  both  given  in  the  laboratory  frame  whereas 

the  elements  of  the  transition  matrix,  are  calculated  in  the  local  frame  and  are  thus  independent  of  the 

orientation  of  the  scattorer;  the  D's,  in  turn,  are  the  irreducible  representations  of  the  rotation  group, i.  e.  the 
matrices  that  transform  the  laboratory  frame  into  the  local  frame.  Equation  (6)  shows  that,  even  when  dealing  with 
an  assembly  of  identical  .scatterers,  we  need  to  calculate  the  transition  matrix  only  once. 

We  have  now  all  the  ingredients  that  are  neces.sary  to  calculate  the  scattered  intensity  from  an  as.sembly  of  identical 
clusters  in  the  vicinity  of  the  reflecting  surface.  It  Is  to  be  stres-sed,  however,  that  the  applicability  of  the  method  of 
the  images  requires  that  we  consider  a  monolayer  of  clusters  that  are  so  arranged  that  the  centers  of  the  corresponding 
spheres  in  different  clusters  have  the  same  distance  from  the  surface.  A  moment’s  thought  will  convince  the  reader 
that  this  condition  implies  that  all  tin'  compound  scatterers  that  are  compo.sed  by  the  actual  objects  and  by  their 
respective  images  are  identical  to  each  other  and  that  the  orientations  of  any  two  of  them  differ  for  a  rotation  around 
an  axis  that  is  perpendicular  to  t  he  surface.  W’ith  this  in  mind  we  can  write  the  total  scattered  intensity  as 
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where  v  and  u'  are  particle  indexes,  N  is  the  number  of  the  particles  on  the  reflecting  surface,  ©)  is  the 

scattered  field  from  a  particle  at  the  position  R  with  orientation  0,  and  the  brackets  denote  the  ensemble  average 
over  both  the  position  and  the  orientation  of  tlie  particles.  The  first  term  on  the  right-hand  side,  the  so  called  self 
term,  contributes  to  the  scattered  intensity  for  any  angle  of  olxservation.  Even  the  second  term,  that  is  a  two-body 
term,  may  become  important  for  nonrandom  distributions  of  the  scatterers,  but,  when  they  are  randomly  distributed 
upon  the  surface,  it  does  contribute  in  the  direction  of  reflection  only.'--  Ultimately,  as  in  this  paper  we  deal  just  with 
monolayers  of  randomly  distributed  scatterers,  we  will  outright  disregard  the  two-body  term  so  that  reduces  to 

7„.  =  yV(/„,)  =  N  J  V(R,0),(.(R,0)dRd0. 


where  t/i(R,  0)  is  the  normalized  density  distribution  function  with  respect  to  the  position  of  the  particles  and  their 
angular  orientation.  Since  we  assume  that  there  is  no  correlation  between  the  position  and  the  orientation  of  the 
scatterers  u>  can  be  factorized  as  w(R,  ©)  =  WR(R)u'e(0);  furthermore  the  distance  of  observation  is  assumed,  as 
usual,  to  be  very  large  with  respect  to  the  size  of  the  reflecting  surface,  so  that  we  can  put  R  =  0  and  write 
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on  account  that 


The  orientation  of  any  two  cIiisIims  iliiri;rs  only  by  a  rotation  around  an  axis  that  is  orthogonal  to  the  surface;  thus,  we 
choose  the  i  axis  of  the  local  reference  frame  that  we  attached  to  each  cluster  to  coincide  with  that  axis,  /?  =  7  =  0, 
and  the  average  involves  the  Euler  angle  o  only.  Now,'" 

=  exp(-j/m 
so  that  eq.  (6)  reads 
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Then,  the  total  scattered  intensity  is 
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where  in  turn 
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ui(a)  =  «>e(«,/?  =  0,7  =  0). 

Equation  (7)  gives  the  scattered  intensity  for  a  general  orientational  distribution  of  the  scatterers,  i.  e.  for  a  general 
t«(a).  In  particular,  if  l(•(o)  =  l/'itr,  i.  e.  when  the  clusters  are  randomly  oriented,  =  6^_o\  this  implies  that  the 
sums  in  eq.  (7)  are  subject  to  the  constraint  ;n  —  in'  —  111"  +  in'"  =  0. 

The  average  over  the  orientations  is,  of  course,  not  necessary  and  is  actually  not  performed  in  two  important 
cases:  when  all  the  compound  .scatterers  are  oriented  alike  or  when  they  have  cylindrical  symtnetry  around  an  axis 
perpendicular  to  the  surface.  The  latter  case  occurs  when  the  centers  of  all  the  spheres  that  form  a  cluster  lie  on  a 
straight  line  that  is  perpendicular  to  the  surface.  However,  even  in  these  two  special  cases  our  orientational  averaging 
procedure  yields  the  correct  results  as  we  will  show  presently.  When  all  the  clusters  are  oriented  alike,  S  can  be 
calculated  with  respect  to  a  local  reference  frame  who.se  ax<»  are  parallel  to  those  of  the  laboratory  frame;  as  a 
consequence  u;(a)  =  ^(0),  where  fi  is  the  Dirac  delta  function,  and  =  1,  so  that  eq.  (7)  becomes 
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If  the  clusters  have  cylindrical  symmetry,  the  elements  not  vanish  for  in  =  m'  only,"*  and  the  same  holds 

true  for  eq.(8),  Itecause  we  just  chose  the  cyliiulrical  axis  as  the  z  axis  of  the  local  frame;  consequently,  eq. 

(7)  simplifies  to 
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on  account  that  2„, =  Iq  =  1  according  to  its  definition.  Thus  even  in  the  above  mentioned  special  cases 
our  orientation  averaging  procedure  works  properly  and  yields  the  trivial  result  that  the  total  scattered  intensity 
equals  .V  times  the  intensity  that  is  scattered  by  an  individual  particle-^'*^. 
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3.  Results  Mild  discussion 


The  theory  tlnit.  we  dc'scrihed  in  (he  preceding  sections  lias  been  tested  against  the  results  reported  by  Johnson 
for  »  single  sphere  near  or  in  contact  with  the  redecting  surface.  Actually,  we  were  able  to  reproduce  with  the 
highest  accuracy  the  rc'snlts  that  Johnson  rcporl.s  in  Figs.  7  -10  of  his  paper.'*  However,  we  are  not  restricted  to 
normal  incidence  so  that  we  are  able  to  calculate  the  full  pattern  of  the  scattered  intensity  at  any  incidence  not 
only  for  a  single  sphere  on  the  reflecting  surface  but  also  for  some  representative  ani.sotropic  model  particles.  In 
fact,  we  perfoi  incd  our  calculations  also  for  a  cluster  of  two  identical  and  mutually  contacting  spheres  lying  on  the 
surface  as  well  as  for  clusters  of  two  identical  muttially  contacting  hemispheres  and  for  linear  chains  of  four  identical 
hemispheres  —  the  neighbouring  hemispheres  arc  in  contact  with  each  other  —  with  their  flat  face  lying  on  the 
reflecting  surface.  With  the  e.xception  of  the  case  of  a  single  sphere,  almost  all  the  patterns  that  we  report  refer  to 
a  dispersion  of  identical  objects  whose  orientations  are  randomly  distributed  as  explained  in  Subsection  2C.  In  our 
calculations  we  assumeil  that  (he  medium  that  fills  the  accessible  half-space  is  the  vacuum  (n  =  1)  and  that  the 
wavelength  of  the  incident  radiation  is  A  =  628.3  nm;  the  refractive  index  of  both  the  spheres  and  the  hemispheres 
is  r»o  =  3  while  the  radius  both  of  the  single  sphere  and  of  the  hemispheres  that  form  a  binary  cluster  was  chosen  to 
be  6,  =  126.0  nm;  in  turn  tiu'  radius  botli  of  the  spheres  that  form  the  binary  cluster  and  of  the  hemispheres  that 
form  the  four-hemispheres  chain  is  />/,  =  100.0  nm.  'fhis  choice  of  the  radii  is  not  casual  as  it  makes  the  total  volume 
of  all  the  clusters,  eiilu-r  of  spheres  or  of  hemispheres,  <‘(|ual  to  the  volume  of  the  single  sphere.  In  fact,  since  we 
are  mainly  interested  in  t  he  effects  of  the  anisotropy,  we  found  it  convenient,  according  to  our  previous  experience,^'* 
to  investigate  objects  with  a  different  geometry  but  containing  the  .same  ipiantity  of  refractive  material.  The  size 
parameter  of  the  single  sphere  (urns  out  to  be  j:,  =  kh,  =  1.26  where.as  the  size  parameter  of  the  spheres  that 
form  the  binary  clusters  is  j'/,  =  F/j/,  =  1.0.  As  a  result,  the  maximum  value  of  /  (  hat  we  had  to  include  into  the 
multipole  expansions,  e<|s.  (1)  and  (3),  was  /,,i/  =  8  for  the  single  sphere,  Lm  =  0  for  (  he  two-hemispheres  clusters 
and  L/if  =  10  both  for  the  t  wo-splu'res  and  the  four-hemispheri's  clusters  in  order  to  achieve  the  convergence  to  four 
significant  digits;  this  convergc'iice  criterion  is.  of  rour.s<‘,  more  strict  that  recpiired  for  an  accurate  display  of  our 
results.  The  need  to  u.se  .so  large  a  value  of  Lm  <?ven  for  the  single  sphere  is  due  to  the  fact  that  our  calculations  do 
not  actually  deal  with  a  single  sphere  but  ratlu'r  with  the  two-.sphere  cluster  that  is  composed  by  the  actual  sphere 
and  by  its  image.  Now,  according  to  Waterman,'^  we  define  the  size  parameter,  x,  of  a  non-spherical  object  as  the 
size  parameter  of  the  smallest  spheri?  that  may  contain  the  object  itself.  Accordingly,  for  the  single  sphere  we  have 
X  =  2.52  and  even  larger  values  for  the  other  clusters  that  we  consider  in  this  paper.  If  one  also  adds  that  the 
refractive  index  is  no  =  3  the  valiK-s  of  Lm  that  we  quoted  above  can  in  no  tvay  be  considered  too  large. 

Before  discussing  our  results  let  us  recall  that  the  theory  of  the  preceding  sections  has  been  set  up  on  the  assumption 
that  the  z  axis  is  orthogonal  to  the  reflecting  surface.  This  choice  is  hardly  compulsory  but,  for  cylindrically-symmetric 
particles,  it  yields  the  automatic  factorization  of  the  matrix  M  as  explained  at  the  end  of  Subsection  2C.  We  also 
stress  that,  according  to  the  set  up  of  the  theory,  the  polar  angles  of  the  direction  of  observation  should  be  in 
the  range  0  <  <ps  <  360®  and  0  <  i?s  <  90®.  Wc  preferred,  instead,  to  display  our  results  with  respect  to  the 
frame  of  reference  that  is  sketched  in  Fig.  1:  the  reflecting  surface  coincides  with  the  xz  plane  and  the  y  axis  is 
thus  orthogonal  to  the  surface  and  directed  towards  the  accessible  halfspace.  Accordingly,  the  scattering  pattern  is 
reported  for  0  <  >ps  <  1^*0®  and  0  <  Os  <  180®.  .Since  the  scattering  pattern  depends  on  the  direction  of  incidence, 
the  latter  need  to  be  specified.  Now,  according  to  eq,  (6),  once  the  matrix  S  is  known,  the  generation  of  the  scattering 
pattern  for  several  values  of  0/  and  is  a  fast  and  low-cost  operation  that  produces  a  large  amount  of  data,  however. 
Therefore  we  resolved  to  report  all  the  scattering  patterns  for  a  single  direction  of  incidence  only.  This  direction, 
that  in  Fig.  1  is  indicated  by  an  arrow,  was  chosen  to  form  an  angle  of  45®  with  the  normal  to  the  surface  and  has 
thus  0i  =  90®  and  (pi  =  225®. 

We  also  chose  to  not  refer  the  state  of  polarization  to  a  pair  of  basis  vectors  that  are  parallel  and  orthogonal  to 
the  scattering  plane,  i.  e.  to  th<?  plane  of  k/  and  ks;  we  preferred  instead  to  project  the  polarization  vector  along 
the  pair  of  unit  vectors  i)  ainl  that  are  tangent  to  the  the  meridians  and  to  the  parallels,  respectively,  of  the  big 
sphere  that  is  depicted  in  Fig.  I.  Of  cours(',  on  account  of  this  choice,  the  appearance  of  cro.s.s-polarization  effects  is 
expected  even  when  the  scatterer  is  a  single  sphere. 

The  quantity  that  we  report  in  Figs.  2-4  as  a  function  of  Os  and  for  the  single  particles  and 

{Ir,q‘) / lo  for  the  dispersions,  where  r  is  (  he  distance  of  observation.  In  the  incident  intensity,  is  the  observed 
intensity,  is  the  orientatioiially-averaged  observed  intensity;  the  indices  i;  and  if  indicate  the  polarization  of 

the  observed  and  of  the  incident  field,  respectively,  and  take  on  the  .symbolic  values  0  and  (p  to  denote  polarization 


along  the  meridians  (I'l-poliiri/alion)  aixl  hIoiik  (lie  panillels  (y7-polari/.alion).  respectively.  Il  should  also  be  noticed 
that  even  when  d.s'  reaches  its  limiting  values,  i)s  =  0®  and  i>s  -  180°,  the  angle  is  still  well  defined  as  this 
angle  characterizes  an  observation  with  a  well  defined  choice  of  the  polarization,  e,  g,  along  the  meridians.  Thus, 
the  limiting  curves  of  our  patterns  {i>s  =  0°  and  i)s  =  180°)  describe  the  observation  of  the  scattered  beani  that 
propagates  along  the  reflecting  surface  at  right  angles  to  the  plane  of  incidence  with  a  polarization  that  depends 
on  (ps-  so,  for  the  i?-polarized  component  of  the  scattered  wave,  when  <pg  =  90°,  is  orthogonal  to  the  surface, 
whereas  when  95  =  0°  or  95  =  180°,  E’'"  is  parallel  to  the  reflecting  surface  and  thus,  as  a  direct  consequence  of  the 
boundary  conditions,  the  scattered  intensity  must  vanish;  for  (he  9-polarized  component  of  the  scattered  wave,  when 
ips  —  0*  or  >ps  =  180®,  E'"'''  is  orthogonal  to  the  surface  whereas,  when  95  =  90°,  E’'^"  is  parallel  to  the  reflecting 
surface  and  thus  the  .scattered  intensity  must  again  vanish.  Since,  for  any  given  polarization,  the  four  extreme  vertices 
of  the  pattern  correspond  to  tin'  same  physical  situation,  the  scattered  intensity  must  have  the  same  value  at  all 
these  extreme  points;  a  further  consequence  is  that  e.  g.  =  0°,9.s  =  0°)  =  =  0°,95  =  90°). 

All  these  fcatnr<*s  can  bi'  seen  in  Kigs.  '2~'i  (  h;i(  report  the  r<-snl(s  of  our  calcidat  ions.  In  particular,  F'ig.  2  reports 
the  scattering  pattern  (.a)  for  a  single  splu're  on  the  surfac<\  (b)  for  a  dispersion  of  randomly  oriented  clusters  of  two 
mutually  contacting  spheres  on  the  surface,  (c)  for  a  dispersion  of  randomly  oriented  linear  chains  of  four  hemispheres 
and  (d)  for  a  linear  chain  of  four  hemispln'ies  along  t  he  z  axis;  the  flat  face  of  all  the  hemispheres  lies  on  the  reflecting 
surface.  The  incident  field  is  9-polarized  and  the  9-polanze<l  component  of  the  .scattered  wave  is  considered.  VVe  first 
remark  that  the  pattern  for  a  single  sphere  in  Fig.  2  (a)  shows  only  a  .symmetry  of  reflection  in  the  plane  of  incidence. 
Therefore,  even  for  such  a  cylindrically  symmetric  object  all  the  directions  of  scattering  must  be  considered  when 
the  incidence  is  not  normal  to  the  reflecting  plane.  In  all  the  four  patterns  in  Fig.  2  we  notice  the  presence  of  a 
strong  scattered  beam  that  propagates  along  the  reflecting  surface  in  the  forward  direction  {0$  =  90°, 95  =  180°). 
We  also  notice  the  two  wings  in  (h<“  pat  tern  from  the  linear  chain  of  four-hemispheres  in  Fig.  2(d).  The  compari.son 
with  the  orientationally  averaged  pati<'rn  in  Fig.  2(c)  suggests  that  the  ob.sorved  wings  are  a  result  of  the  particular 
orientation  of  the  scatterers.  bet.  ns  now  consider  (In'  limiting  curves  at  tig  =  0°  and  Og  =  180°.  It  is  quite  evident 
that  these  curves  are  almost  flat  for  all  the  clusters  but  are  in  no  way  flat  for  the  single  sphere.  This  suggests  that 
the  observation  of  the  scattered  light  that  propagatis;  along  the  reflecting  surface  may  yield  useful  information  on 
the  shape  of  the  scattering  particles. 

This  is  confirmed  by  (he  pa(t<'rns  in  Fig.  d  that  report  our  results  (a)  for  a  single  sphere  and  for  the  a.ssemblies 
(b)  of  randomly  oriented  two-spln'i-os  clusters,  (c)  two-hemi.spheros  clusters  and  (d)  four-hemispheres  linear  chains. 
The  incident  wave  is  d-polarized  ;md  the  9-polarized  component  of  the  .scattered  wave  is  considered,  .so  that  all  the 
patterns  describe  cross-polarization  effects,  bet  us  remark  fir.st  that  the  patterns  in  (a)  and  (b),  that  refer  to  a  single 
sphere  and  to  the  two-sphere  clusters,  pro.sent  a  couple  of  peaks  for  non-limiting  values  of  ^g  and  95:  such  peaks  are 
not  present  in  the  patterns  in  (c)  and  (d),  that  refer  to  the  two-hemispheres  and  to  the  four-hemispheres  clusters. 
Furthermore,  the  patterns  for  all  the  clusters,  either  of  spheres  and  of  hemispheres  reach  titeir  maximum  value  at 
the  four  vertices  (ii)g  =  0°,  180°  and  pg  =  0°,  180°);  this  is  not  the  case  for  the  single  sphere  in  (a).  The  differences 
of  the  patterns  of  the  single  sphere  and  of  tin'  clusters  can  be  attributed  on  one  hand  to  the  cylindrical  symmetry 
of  the  single  sphere  and  on  the  other  hand  to  (he  anisotropy  of  the  cln.stcrs  in  spile  of  their  random  orientation.  In 
other  words,  the  process  of  averaging  does  not  cancel  the  effects  that  are  due  to  the  ani.sotropy  of  the  scatterers.  The 
patterns  in  Fig.  3  suggest  that  the  effects  of  the  anisotropy  are  still  visible  provided  that  the  ob.servation  is  made 
with  the  appropriate  polarization. 

In  Fig.  4  we  report  the  patterns  for  the  same  objects  that  we  considered  in  Fig.  3;  the  only  difference  is  the  choice 
of  the  polarization,  i.  e.  the  incident  wave  is  9-polarized  and  the  iJ-polarized  component  of  the  scattered  wave  is 
considered.  The  patterns  in  (a)  and  (c),  that  refer  to  the  single  sphere  and  to  the  two-hemispheres  clusters,  show 
a  strong  similarity;  such  a  similarity  is  also  shown  by  the  patterns  in  (b)  and  (d)  that  refer  to  the  two-spheres  and 
to  the  four-hemispheres,  respectively.  Thus,  the  general  structure  of  the  patterns  seems  to  be  related  to  the  number 
of  the  spheres  that  compose  the  actual  scatterers,  providerl  that  two  hemi.sphercs  are  counted  a.s  one  sphere.  In  this 
respect  we  recall  that  all  the  scatterers  that  we  investigate  in  this  paper  have  the  same  volume. 

We  do  not  report  any  result  for  (  he  case  in  which  both  the  incident  and  the  .scattered  wave  are  considered  to  be 
«>-polarized;  these  patterns  do  not  add  in  fact  any  information  worth  of  a  separate  comment. 


4.  Conclusions 


The  results  that  wo  ropoi  lod  in  l  lio  Section  ;{  siig^fost  a  luimhor  of  considerations  both  on  the  applicability  of  the 
small  particle  approximation  to  the  problem  at  hand  and  to  the  ability  of  onr  approach  to  yield  useful  information 
for  the  interpretation  of  the  experimental  data. 

In  the  present  problem  the  range  of  applicability  of  the  Rayleigh  scattering  approximation  (RSA)  is  even  more 
restricted  than  in  the  ca.se  of  the  isolated  particles.  The  n.se  of  the  theory  of  the  images,  indeed,  implies  that  one 
has  to  deal  with  the  compound  .sentterer  that  is  composed  by  the  actual  particle  and  by  its  mirror  image.  According 
to  the  definition  that  we  (pioted  in  the  pn-ceding  .section,*®  the  effective  size  parameter  of  such  a  compound  object 
is  larger  than  the  size  parameter  of  tlie  act  ual  scalterer.  Of  course,  if  the  actual  particle  is  not  in  contact  with  the 
reflecting  surface  the  ('fh  etivi'  size  paraiiK'ier  is  bound  to  increa.se.  Obviously,  the  RSA  is  in  no  way  applicable  to 
the  scalterers  that  we  consider  in  l  itis  paper  bot  h  because  the  size  parameter  of  each  component  sphere  (and  of  each 
hemisphere  too)  is  not  small  an<l  becaii.si-  the  ntfractive  ind<\\  il..s<,df  is  not  close  to  unity.  In  this  respect  we  recall 
that  we  had  to  consider  mnll  ipole  fields  up  to  l,\i  =  It),  in  contrast  to  the  RSA  that  amounts  to  consider  the  dipole 
field  only,  i.  e.  /,.»  =  I 

The  approach  that  wt'  n.sc'd  in  this  paper  has  been  designe<l  just  to  deal  with  the  scatli.'ring  properties  of  anisotropic 
particles  without  resorting  to  any  approximation  except  for  the  truncation  of  the  mnltipole  expansions  whose  effect 
can  easily  be  checked,  however.  The  nsefulin'ss  of  onr  approach  stems  from  il.s  goneraSty  that,  in  turn,  derives 
from  the  transformation  properties  under  rotation  of  the  spherical  mnltipole  fields  and  of  the  transition  matrix.  As 
a  result,  we  were  able  to  calculate  with  a  limited  computational  effort  th<'  scattered  intensity  from  assemblies  of 
randomly  oriented  anisotropic  scatterers  on  the  reflecting  surface. 

The  results  of  our  calculations,  in  particular  the  n'sidts  in  Figs.  2  and  3,  show  that  the  .scattering  pattern  of  onr 
model  anisotropic  scatterers  an-  remarkably  different  from  those  of  the  spherical  particles.  Therefore,  provided  this 
behavior  pertains  also  to  more  general  atii.sotropic  scatterers  than  the  model  .scatterers  that  we  considered  so  far,  a 
polarization  analysis  of  the  grazing  scatt<Med  ititensil.y  could  yield  n.sefni  information  on  the  possible  anisotropy  of 
the  particles  on  the  surface. 


Appe.»dix  A:  Scattering  from  a  cluster  of  spheres 

We  consider  a  cluster  of  N  spheres  of  radius  ha  and  (po.ssibly  complex)  refractive  index  tIq  whose  centers  lie  at  R,,; 
the  cluster  is  assumed  to  be  embedded  into  a  homogeneous  medium  of  refractive  index  n.  The  incident  field,  that  we 
assume  to  be  the  plane  wave  whose  mnltipole  expansion  is  given  by  eq.  (1),  excites  within  each  sphere  the  internal 
field 

plm 

where  r®  =  r  —  Ro-  By  writing  the  scattered  field  as  a  superposition  of  the  fields  that  are  scattered  by  the  single 
spheres, 

o  plm 


the  unknown  amplitudes  can  be  calculated  by  imposing  the  boundary  conditions  acro.ss  the  surface  of  each 

one  of  the  spheres  in  the  cluster.  In  this  respect  we  notice  that  the  scattered  field  is  written  in  terms  of  multipole 
fields  that  are  referred  to  different  origins  wherea.s  the  incident  field  is  written  in  terms  of  mnitipoles  that  are  referred 
to  the  origin  of  the  frame  of  reference.  Therefore,  before  the  boundary  conditions  are  applied  we  have  to  use  an 
appropriate  addition  theorem  to  refer  all  the  mnltipole  fields  to  the  same  origin.*®  For  the  incident  plane  wave,  e.  g., 
this  theorem  yields 

Im 


where 


The  quantities  '<i‘<'  iIk'  i  Ii'iik'hIs  of  tlio  niatiix  ihal  peifoiiiis  ilio  iiansrer  of  origin  of  tin;  miill.ipolo  (icids 

from  the  unique  origin  ni  Ro  lo  R„;  of  lotirso,  in  lln>  present  case  the  unique  origin  Ro  coincides  with  the  origin  of 
the  reference  frame,  so  I  hat  Rd  =  0.  VVe  hav<> 


7(ry)  _  i  ,  _  ; 


^I'v'  -  '  ^  E  ^  ~ 

J 


x(7/+i  ( A',  R„  nK'(l. w'  +  //). 

where  R^o  =  Rn  -  R.i  ''iid 

G, A'.  R.)  =:  ■lT5];;'-''^"/,{/'.m';/,m)j,(A/0V;,„._,.,.(R).  (A2) 

A 

In  the  preceding  tapial  ions  the  (  "s  are  the  ( 'lel>sch-(<ordan  coellicients  and  the  quantities  lx  are  the  (laiint  integrals'  * 
that  are  (lefined  as 


Ix(/',m‘;Am)  -  J 


_  <-•(/,  A, /^O,  OK  •(/,  A, -  m). 

V  I^Ci/'  +1) 

We  will  also  mvd  the  inverse'  transfer  from  0,,  to  Rn  in  orth'r  to  get  t  he  scat  tere'd  Reid  in  terms  of  mnitipoles  referred 
to  an  unique  origin,  I'his  transfi'r  is  effected  hy  a  matrix  whose  elemi'nts,  ,  are  idimtical  to  the  elements 

*^a(m  I'm'  sl’ove  except  t  hat  in  (!  oin'  must  .suhstitnti-  Ro,,  for  R„o, 

Anyway,  once  the  incident  and  the  .scatten-d  field  are  expressed  in  terms  of  multipoles  referred  in  turn  to  eacli  one 
of  the  sites  of  the  cluster  and  the  hoiindary  conditions  arc  imposed,  one  gets,  for  each  p,  o,  /  and  m,  four  equations 
among  which  the  internal  amplitudes  cun  be  eliminated.  As  a  result,  the  scattered  amplitudes  turn  out 

to  be  the  solution  of  the  system  of  linear  non- homogeneous  equations 

MA  =  -W,  (A3) 

where  Ad  =  7v~'  -»->f  and  A  and  >V  are  tin'  one-coliimn  vectors  of  the  amplitudes  '^^Tim'  rcs|>'5clively.  V. 

is  a  diagonal  matrix  of  elements 

^o*/m,o' I'm'  “  Aio'/fo/ 


„  (p)  _  (1  +  n„i>,,t)iii(v„kb„)ii',{vhh„)  -  ( I  -f  >in^.;)n{(»o^fto)«((»^’o) 

(1  -I- rin/ipi )i'/(»a*d>„)ii)[(rU-6„)  -  (I  -I- ii„/i,,-.)«|(naA6o)ii>((iiA6„) ’ 


ff„  =  —  -  I,  i//(r)  =  J-jiij-).  Wi(j)  =  .r/ij"(r). 

Ra\'^  and  fJoJ’’  coincirle  with  the  Mie  coelficienls  hi  and  <i/,  re.spectively,  for  the  scattering  hy  a  homogeneous  sphere 
of  radius  6^  and  refractive  index  ii„  emheddi'd  in  a  homogeneous  medium  of  refractive  index  ri,  The  last  matrix  that 
we  need  to  define  is  W  whosi'  elements  are 
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X 


=  (I  -  V  -  (I 


^C(l,/  +  I  -  (<,,,, '.I-,  -/(.  Ill  +  /<  •'/+ 1  -  A,, .  n.,1  I./'.  III'  +  /<). 

h 

where  G/m,J'm'(A  .  R-o  <.')  is  still  given  hy  eci.  (A‘2)  exeopl  for  the  siihstitulion  of  /i\'*  to  and  R.ao'  =  Ro'  -  R-a- 
The  formal  solution  of  <'(|.  (AM)  is 

A  =  -Af-'VV 


This  equation  is  rather  similar  to  eq.  (1)  lull  vvi'  are  |>r<‘V'enre<l  to  identify  Ad”'  with  S  hy  the  fart  that  both  the  >t’s 
and  the  W’s  are  tlu'  amplitudes  of  imiltipoles  that  are  referred  to  different  origins.  However,  multiplication  of  the 
preceding  equation  hy  and  summation  over  o,  />,  I  and  in,  after  substitution  for  the  W’s  of  their  expression 

in  eq.  (Al)  yield.* 

aU')  _  V  V  7'''''''  A*'"’' 

^Itn  ~  /  .  /  -  I'm' ' 

and,  ultimately, 


ctf'P')  _  V' V'  V'  f  M-'l*''''"*  -/U'.p'i 

^Im.l'm'  ~  2-.I  ''l"i  '>  1^1  I'  '  In  I.M.,,’ I  I  'M'.l'i"'' 


<1 '  ifl.  M  n*  !.*  M' 
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Fig.  I.  SkcU'll  of  till'  jn we  <  Iiom'  for  llir  ilis|>l:iv  of  oin  ii'siih>.  I'lic  tcflci  I iiig  Miifarc  t oiiicitlcN  »vi(li  llir  xz 
ftxin  and  the  y  axis  |>oiiils  lowanK  l  lic  atx'ssililc  half  s|>ai  r  who.si'  u  liat'livr  imlcx  is  »  =  I .  I'lu'  jxilar  angles  of  llie  direelioii 
of  observation  range  in  iln-  inleivals  It”  <  iK  <  |sii“  ami  il“  <  y- .  <  |S(l°  The  <lin'i  lion  of  in< idem c,  fe;,  llial  is  imih  aleil 
by  the  arrow,  has  llu'  (xtlai  angles  il/  ~  'ltl°  ami  y-/  —  I'lx-  airangeineiil  of  the  four. hemispheres  linear  chain  whose 

scattering  pattern  is  repot  leil  in  I'ig  Jftl)  below  is  also  shown 


Fig.  2.  Pattern  of  the  scallered  inlensily  (a)  for  a  single  sphere  on  the  rrdlecting  surface,  (b)  for  a  dispersion  of  ran¬ 
domly-oriented  two-spheres  clusters,  (c)  for  a  dispersion  of  randoinlv  oriented  fonr-heinisphcres  linear  chains  and  (d)  for  a 
four-hemispheres  linear  chain  oriented  as  shown  in  Pig.  I.  The  wavelength  of  the  incident  wave  is  A  =  028. 3  nm  and  the 
refractive  index  of  all  the  scalterers  is  iiti  =  3.  The  radins  of  the  single  sphere  in  (a)  is  ft,  =:  l2C.Uniti  rvhereas  the  radius  both 
of  the  spheres  of  the  two-spheres  cinsler  in  (b)  and  of  the  hemispheres  of  the  four-hemispheres  linear  chain  in  (c)  ami  (d)  is 
it,  =  lOOntn.  We  artnallv  report  (in  stpiare  meters)  t he  ipiantities  r'l^.^./lu  in  (a)  and  (d)  and  in  (b)  and  (c)  as 

a  function  of  the  angh's  of  observation  ds  and  yss.  is  the  rlistanx-  of  observation,  /q  is  the  incident  intensity,  /y..  is  the 
observed  intensity  and  (/,  ,.}  the  orient.itionallv-averaged  ob.serxxd  intensity. 

Fig.  3.  Pattern  of  the  s(attereil  intensity  (,i)  for  a  singh'  sphere,  (b)  for  a  dispersion  of  randomly  oriented  two-spheres 
clusters,  (c)  for  a  ilispersion  of  ramlomly-oriented  two-hemisplx  res  <  bisters  ami  (d)  for  a  dispersion  of  randomly-oriented 
font-hemispheres  linear  ih.iins  I'he  radii  .iml  the  refractive  imixes  are  as  in  I'ig.  2:  the  radius  of  lh<^  hemispheres  of  the 
two-hemispheres  cluster  ei|ii.ds  the  r.idiiis  of  tlx-  single  sphere.  We  report  (in  .stpiare  meters)  l^.fl/la  in  (a)  and  r^(/y.a)//o 
in  (b),  (c)  and  (d). 


Fig.  4.  Same  as  I'ig  1  except  that  i^/,i,.//n  and  r^{l,i,)/la  are  consideretl. 
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